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Abstract 

The global data on solar neutrino rates and spectrum, including the SNO 
charged current rate, can be explained by LMA, LOW or the energy indepen- 
dent solution - corresponding to near-maximal mixing. All the three favour 
a mild upward renormalisation of the CI rate. A mild downward shift of the 
B neutrino flux is favoured by the energy independent and to a lesser extent 
the LOW solution, but not by LMA. Comparison with the ratio of SK elastic 
and SNO charged current scattering rates favours the LMA over the other 
two solutions, but by no more than 1.5cr. 



The suppression of the solar neutrino flux has been confirmed now by a 
number of experiments [1-4], covering a neutrino energy range of 0.2 — 20 
MeV. However there is still considerable uncertainty regarding its energy de- 
pendence. An energy independent suppression factor has been advocated by 
several groups over the years [5]. This was shown to be disfavoured how- 
ever by the combined data on the total suppression rates [6]. More recently 
the weight of experimental evidence changed in favour of an energy inde- 
pendent solution following the SK data on the day/night spectrum, showing 
practically no energy dependence nor any perceptible day/night effect [3]. 

It was recently shown by us in [7] that an energy independent solution can 
describe the global data on rates and spectrum satisfactorily, with reasonable 
readjustments of the CI experiment rate and the neutrino flux. This has 
been corroborated now by other groups [8,9]. It was also suggested in [7] 
that the measurement of the charged to neutral current ratio by the SNO 
experiment will be able to distinguish the energy independent solution from 
the conventional MSW solutions to the rates and spectrum data, i.e. the so 
called LMA and LOW solutions. Recently the SNO experiment has produced 
its first data [4] , providing the charged current scattering rate with reasonable 
precission. In the absence of the neutral current data one can combine the 
CC rate of SNO with the elastic scattering rate of SK to distinguish the 
above solutions from one another. The purpose of the present exercise is two 
fold. Firstly we reevaluate the experimental status of the energy independent 
solution vis a vis the LMA and LOW solutions after the inclusion of the 
SNO data. Secondly we plug an important gap in the analysis of [7], which 
arose out of neglecting the coherent term in the solar neutrino oscillation. 
Inchision of this term allows us to probe the parameter space down to Am^ ~ 
10-11 eV^ i.e. the vacuum oscillation region, which was not possible in [7]. 
We also take this opportunity to update the result using the full SK data 
sample of 1258 days, which became available more recently [3]. We shall see 
that the conclusion of [7] are not changed by this and the inclusion of the 
SNO CC rate in the global rates and spectrum data. On the other hand a 
direct comparison with the ratio of SNO (CC) and SK (ES) scattering rates 
does favour the LMA solution over the energy independent (and to a lesser 
extent the LOW) solution, but at only ^ 1.5a level. 

Table 1 shows the suppression rate or survival probability of the solar 
neutrino (P^,^,) from the combined Ga [1], CI [2], SK [3] and SNO [4] exper- 
iments along with their threshold energies. The corresponding compositions 
of the solar neutrino flux are also indicated. The SK survival rate shown in 
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parantheses is appropriate for the oscillation of into another active flavour 
('^^.r)! on which wc shall concentrate here. It is obtained by subtracting the 
neutral current contribution of v^^r from the SK rate. All the survival rates 
are shown relative to the standard solar model (SSM) prediction of BPBOO 



The apparent energy dependence in the survival rates of Table 1 is con- 
ventionally explained in terms of the vacuum oscillation (VO), small and 
large angle MSW (SMA and LMA) as well as the LOW solutions [11]. The 
VO and SMA solutions, showing strong and nonmonotonic energy depen- 
dence, are essentially ruled out now by the SK day/night spectrum [3]. On 
the other hand the spectrum data is compatible with the LMA and LOW 
solutions, which predict modest and monotonic decrease of the survival rate 
with energy. 

We shall fit the combined data on the survival rates and the SK spectrum 
with the energy independent solution as well as the four traditional solutions 
mentioned above. In order to reconcile the energy independence of the spec- 
trum with the apparent energy dependence in the rates of Table 1, we shall 
consider the following variations in the CI rate and B neutrino flux, since the 
CI experiment [2] has not been calibrated while the B neutrino flux is very 
sensitive to the underlying solar model. 

i) We shall consider an upward renormalisation of the CI rate by 20% (i.e. 
2(t), which will push it marginally above the SK and SNO rates. This 
is favoured not only by the energy independent solution but also the 
LMA and LOW solutions, showing monotonic energy dependence. 

ii) We shall consider a downward variation of the B neutrino flux of BPBOO 



by upto 2a. A downward renormalisation of this flux is favoured by the 
energy independent solution and to a lesser extent by LOW, though not 
by the LMA solution. It is also favoured by some helioseismic models, 
e.g. the model of [12] giving /b = (4.16 ± 0.76) x 10^/cm^/sec. 

The deflnition of used in our flts is 
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where i, j run over the number of data points, and both the experimental and 
theoretical values of the fitted quantities are normalised relative to the BPBOO 
[10] prediction. The experimental values for the total rates are the ones 
shown in Table 1, while the SK day/night spectra are taken from [3]. The 
error matrix atj contains the experimental errors as well as the theoretical 
errors of the BPBOO predictions along with their correlations. The latter 
includes the uncertainty in the B neutrino flux of eq. (|l|). The error matrix 
is evaluated using the procedure of [13]. The details of the solar code used 
is described in [14]. As in [7] we vary the normalisation of the SK spectrum 
as a free parameter to avoid double counting with the SK data on total rate. 
Thus there are (2 x 19 — 1) independent data points from the SK day/night 
spectrum along with the 4 total rates giving a total of 41 points. In addition 
to the best-fit parameter values and xLin 'we shall present the goodness of 
fit (g.o.f.) of a solution, which represents the probability of a correct model 
having a ^ this Xmm- Finally for the oscillation solutions with matter 
effects we shall also delineate the 90%, 95%, 99% and 99.73% (3a) allowed 
regions in the mass and mixing parameter (Am^ — tan^ 6) plane. These 
regions correspond to ^ Xmm + ^X^y where Ax^ = 4.61, 5.99, 9.21 and 
11.83 respectively for two parameters and Xmin is the global minimum. 

Energy Independent Solution: 

Table 2 summarises the results of fitting the global rate + spectrum data 
with an energy independent survival probability 

Pu.u. = l-\sm'2e. (3) 

It shows that even without any readjustment to the CI rate or the B neutrino 
flux the energy independent solution gives an acceptable g.o.f. of 24%. An 
upward renormalisation of the CI rate by 20% improves the g.o.f. to 42%. 
And varying the B neutrino flux downwards improves it further to 49%, 
which corresponds to a renormalisation factor Xb = 0.7 for the B neutrino 
flux. Note however that the g.o.f. and the best-flt value of the mixing angle 
for Xb = 1 are very close to those for Xb = 0.7. This is because the large 
error in the B neutrino flux of eq. (|I|) is already incorporated into the error 
matrix of eq. (0). It automatically chooses a lower value of flux, close to 
Xb = 0.7, even without floating this parameter. 

Traditionally the energy independent solution is associated with the 
averaged vacuum oscillation probability at distances much larger than the 
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oscillation wave-length, as originally suggested by Gribov and Pontecorvo 
[15]. As we shall see below however an effectively energy independent solution 
holds around the maximal mixing region over a wide range of Am? even after 
including all the matter effects. 

Region of Energy Independent Solution in the Am^ — tan^ 9 Plane: 

The energy dependence of the survival probability arises from different 
sources in different regions of the parameter space. 

i) For Am? / E < 10"^'' eV the earth regeneration effect can be safely 
neglected. Then the survival probability 

Pv.v. = Pi cos' e + P2 sin' 6 + 2/p^sin 9 cos (^^^ j , (4) 

where L is the distance between the neutrino production point at the 
solar core and its detection point on earth; and P2 (= 1 — Pi) is the 
probability of the produced v^. emerging from the sun as the heavier 
mass eigen-state 

z/2 = z^e sin6' + z/^ cos6'. (5) 

The coherent interference term, represented by the last term of eq. (^, 
is responsible for the nonmonotonic energy dependence of the socalled 
Just-So (VO) solution. 

ii) For Am? jE ~ 10^^^ — 10^^^ eV, the coherent term is negligible while 
the earth regeneration contribution can be significant. Besides over 
a large part of this region, represented by the MSW triangle, the v^, 
is adiabatically converted into V2 in the sun, i.e. P2 = 1. Thus the 
day/night averaged probability 

= ^ + 4(1 -2,, cos 2^ + , I)' 

where 

Here p is the matter density in the earth and Ye the average number 
of electrons per nucleon [16]. The regeneration contribution is always 
positive and peaks around r]E ~ 1, i-e. Am? jE ~ 3 x 10^^^ eV. This 
is responsible for the LOW solution. 
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iii) Finally for ^m? /E > 10~^^ eV the survival probability can be approx- 
imated by the average vacuum oscillation probability of eq. @. The 
MSW solutions (LMA and SMA) lie on the boundary of the regions ii) 
and iii), i.e. Am^ ~ 10~^ eV^ for E 1 MeV. The survival probability 
PueUe goes down from 1 — | sin^ 29 (> 0.5) to sin^ 9 (< 0.5) in going up 
from Ga to SK (SNO) energy. 

All the solar neutrino rates except that of SK have ^ 10% error, which is 
also true for the SK energy spectrum. The SK normalisation has at least sim- 
ilar uncertainty from the B neutrino flux. Therefore we shall treat solutions, 
which predict survival probability Pu^u^^ within 10% of eq. over the Ga to 
SK energy range, as effectively energy independent solutions [7]. Moreover 
the predicted Ga, CI and SNO rates will be averaged over the respective 
energy spectra, while the predicted SK rates will be averaged over 0.5 MeV 
bins, corresponding to the SK spectral data, since experimental information 
is available for only these averaged quantities. 

Fig. 1 shows the region of effectively energy independent solution as per 
the above definition. The parameter space has been restricted to Am^ < 
10~^ eV^ in view of the constraint from the CHOOZ experiment [17]. One 
sees that the energy independent solution (^) is effectively valid over the 
two quasi vacuum oscillation regions lying above and below the MSW range. 
Moreover it is valid over a much larger range of Am^ around the maximal 
mixing region, since the solar matter effect does not affect Pu^ue tan^ ^ = 1. 
It is this near-maximal mixing strip that is relevant for the observed survival 
probability, Pu^ue ~ 1/2. The upper strip (Am? = 10^^ — 10~^ eV^) spans 
the regions iii) and part of ii) till it is cut off by the earth regeneration effect. 
The lower strip (Am^ = 10^^ — 5 x lO^^*' eV^) spans parts of region ii) and 
i) till it is cut off by the coherent term contribution. Note that this near- 
maximal mixing strip represents a very important region of the parameter 
space, which is favoured by the socalled bimaximal mixing models of solar 
and atmospheric neutrino oscillations [18,19]. 

One can easily check that averaging over the SK energy bins of 0.5 MeV 
has the effect of washing out the coherent term contribution for Am^ > 
2 X 10~^ eV^. Hence the contour of effectively energy independent solution, 
shown in ref. [7] neglecting the coherent term, is correct. But including this 
term enables us now to trace the contour down to its lower limit. It was 
claimed in ref. [9] that the lower strip disappears when one includes the 
coherent term contribution. This may be due to the fact that their predicted 
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rate in the SK energy range was not integrated over the corresponding bin 
widths of 0.5 MeV. 

To get further insight into the oscillation phenomenon in the maximal 
mixing region we have plotted in Fig. 2 the predicted survival rates at max- 
imal mixing against Am^ for the Ga, CI, SK and SNO experiments. In 
each case the rate has been averaged over the corresponding energy spec- 
trum. This is similar to the Fig. 7 of Gonzalez-Garcia, Pena-Garay, Nir 
and Smirnov [16]. As in [16] the predictions have been shown relative to the 
central value of the B neutrino flux of BPBOO along with those differing by 
±20% from the central value, i.e. Xb = 1 ± 0.2. We have found that these 
curves are in good agreement with the corresponding ones of [16]. The two 
regions of < 10% energy dependence are indicated by vertical lines. One 
can easily check that in these regions the central curves lie within 10% of the 
energy independent prediction R = 0.5 (only the SK rate is higher due to the 
neutral current contribution). One can clearly see the violent oscillation in 
the Just-So (VO) region at the left and the LOW solution corresponding to 
the earth regeneration peak of the Ga experiment in the middle. The LMA 
and SMA are not identifiable since they do not occur at maximal mixing 
angle. It should be noted that the gap between the two energy independent 
regions due to the earth regeneration effect in Fig. 1 is a little narrower than 
here. This is due a cancellation between the positive contribution from the 
earth regeneration effect and the negative contribution from the solar matter 
effect at tan^ ^ < 1. It ensures that the resulting survival rate agrees with 
the energy independent solution over a somewhat larger range of Am^ 
slightly below the maximal mixing angle. 

The observed rates from the Ga, CI, SK and SNO experiments are shown 
in Fig. 2 as horizontal lines along with their la errors. With 20% downward 
renormalisation of the B neutrino flux {Xb = 0.8) the energy independent 
prediction is seen to agree with the SK rate and also approximately with 
SNO. It is higher than the CI rate by about 2a. The agreement with the 
Ga rate can be improved by going to a little smaller 6 and compensating 
the resulting deviation from the other rates by a somewhat smaller Xb as 
in Table 2. Nonetheless the maximal mixing solution for Xb = 0.8, shown 
here, is in reasonable agreement with the observed rates over the energy 
independent regions. The earth regeneration effect can be seen to improve 
the agreement with the Ga experiment for the LOW solution. 

The SMA, LMA, LOW and VO Solutions: 
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Tables 3 and 4 summarises the results of fits to the global rates + spec- 
trum data in terms of the conventional oscillation solutions. Table 3 shows 
solutions to the data with observed and renormalised CI rate with the neu- 
trino flux of BPBOO {Xb = 1), while Table 4 shows the effects of renormalis- 
ing this B neutrino flux downwards by 25% {Xb — 0.75). The corresponding 
90%, 95%, 99% and 99.73% (3(7) contours are shown in Fig. 3. 

As we see from these tables and Fig. 3 the SMA solution gives rather 
poor fit in each case, with no allowed contour at 3a level. This result agrees 
with the recent fits of [20] to the global data including SNO. The fit of [21] to 
these data shows a small allowed region for SMA solution at 3a level due to 
a slightly different method of treating the normalisation in the SK spectrum 
data, as explained there. 

The LMA and LOW solutions give good fits to the original data set, 
which improve further with the upward renormalisation of the CI rate by 
2a. This is because the monotonic decrease of rate with energy, implied by 
these solutions, favours the CI rate to be marginally higher than the SNO 
and SK rates as mentioned earlier. For the renormalised CI case, downward 
renormalisation of the B neutrino fiux by 25% is seen to give a modest 
increase (decrease) of g.o.f. for the best LOW (LMA) solution. On the 
other hand the allowed ranges increase in both cases as we see from Fig. 3b 
and c. Combined together they imply that the downward renormalisation 
of the B neutrino fiux modestly favours the LOW solution but makes little 
difference to the LMA. Its main effect on these two solutions is increasing 
their allowed ranges in the parameter space. Comparing Fig. 3c with Fig. 1 
shows that much of these enlarged ranges of validity correspond to effectively 
energy independent solution. It is intersting to note that the best-fit values of 
parameters in the LMA region are same for CI observed and CI renormalised 
cases while the Xmin improves for the latter. This shows that the best-fit 
already chose a probability at CI energy, which is a little higher than that 
at SK/SNO energy. Renormahsing the CI rate brought that point up to the 
fitted curve without changing the best-fit parameters. 

While the best vacuum solution seems to show remarkably high g.o.f. 
particularly for renormalised CI rate and B neutrino fiux, its regions of va- 
lidity are two miniscule islands just below the lower energy independent strip 
of Fig. 3b, c. This solution has also been obtained in the global fits of ref. 
[20,21] as well as the SK fit to their rate and spectrum data [3]. The position 
and size of its range of validity suggest this to be a downward fiuct nation 
of the effectively energy independent quasi vacuum oscillation rather than a 



8 



genuine VO solution of Just-So type. To get further insight into this solution 
we have analysed the resulting energy dependence. It shows practically no 
energy dependence below 5 MeV, but a 15% fall over the 5 — 12 MeV range. 
The latter seems to follow the SK spectral points rather closely amidst large 
fluctuation. To check the stabihty of this trend we have repeated the flt 
to the SK spectral data points, plotted over 8 broad energy bins shown in 
[3], which show much less fluctuation than the 2 x 19 points sample. The 
solution completely disappears from this fit. This confirms that the above 
VO solution is simply an artifact of the sampling of the SK spectral data. 

For completeness we summarise in Table 5 the best fits of the above 
solutions with free B neutrino flux normahsation. The SMA solution favours 
a very low B neutrino flux {Xb — 0.5), which raises the SK and SNO rates 
more than the CI, thus accentuating the nonmonotonic energy dependence 
of Table 1. Still the g.o.f. of the SMA solution is rather marginal. On the 
other hand the LMA solution favours Xb > 1, which suppresses the SK and 
SNO rates more than the CI, resulting in a monotonic decrease of rate with 
energy. But the corresponding g.o.f. are no better than those of Table 3. 
The results of the LOW and VO fits are similar to those of Table 4. 

Compcirison with the Ratio of SK and SNO Rates: 

With the measurement of both charged and neutral current scattering 
rates at SNO it will be possible to discriminate between the above solutions, 
since the B neutrino fiux factors out from their ratio [7]. In the absence of 
the neutral current data from SNO one can try to make a similar comparison 
with the ratio of SK elastic and SNO charged current scattering rates, 

R'Ij, = XbPu^u^ + r(l - P.^.JXb, r = a^Jacc ^ 0.17, (8) 

R'sNO ^ ^bPu^i^^, (9) 

where we have assumed a common survival rate neglecting the small differ- 
ence between the SK and SNO energy spectra [22]. One can eliminate Pu^i,^ 
from the two rates; and the resulting B neutrino fiux can be seen to be in 
good agreement with the BPBOO estimate [4]. Alternatively one can factor 
out the flux from the ratio 

RsK/RsNO = '^-r + r/P.^.^. (10) 

Table 6 shows the best fit values of the above ratio for the LMA, LOW and 
the energy independent solutions along with the corresponding predictions 



9 



for R'sNo/ -^^'no ~ PfeVe- '^^^ predictions of the maximal mixing solution 
is also shown for comparison. While the predictions for the CC/NC ratio 
differ by ~ 50% they differ by only about ~ 15% in the case of the ES/CC 
ratio. The observed ratio R'sk/^'sno ^^^^ favour the LMA over the 
LOW and energy independent solutions; but even the largest discrepancy is 
only ~ 1.5cr. With the expected sample of several thousand CC and NC 
events from SNO one expects to reduce the la error for each of these ratios to 
about 5%. Then one will be able to discriminate between the three solutions 
meaningfully, particularly with the help of the CC/NC ratio from SNO. 

S.C. and S.G. would like to acknowledge Abhijit Bandyopadhyay for dis- 
cussions. 
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Table 1: The ratio of the observed solar neutrino rates to the corresponding 
BPBOO SSM predictions. 

"experiment BPmo composition 

CI 0.335 ± 0.029 B (75%), Be (15%) 

Ga 0.584 ± 0.039 pp (55%), Be (25%), B (10%) 

SK 0.459 ± 0.017 B (100%) 
(0.351 ±0.017) 

SNO(CC) 0.347 ± 0.027 B (100%) 



Table 2: The best-fit value of the parameter, the x^m the g.o.f from a 
combined analysis of rate and spectrum with the energy independent solution 
given by eq. (3). 









/ tan^ 9 \ 








Xb 


sin^ 29 


or 


A mm 


g.o.f 








V cot^ 9 J 






Chlorine 


1.0 


0.93(0.57) 


46.06 


23.58% 


Observed 


0.72 


0.94(0.60) 


44.86 


27.54% 


Chlorine 


1.0 


0.87(0.47) 


41.19 


41.83% 


Renormalised 


0.70 


0.88(0.48) 


38.63 


48.66% 
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Table 3: The best-fit values of the parameters, the Xiuin ^^^^ ^^'^ g-O-f from 
a combined analysis of the CI, Ga, SK and SNO CC rates and the SK day- 
night spectrum in terms of Ve oscillation into an active neutrino, including 
the matter effects. Xb is kept fixed at the SSM value (=1.0). 





Nature of 
Solution 


in eV^ 




tan^^ 


■ 


g.o.f 




SMA 


5.28 X 10" 


-fi 


3.75 X 10"^ 


51.14 


9.22% 


CI 


LMA 


4.70 X 10" 


-5 


0.38 


33.42 


72.18% 


Obsvd. 


LOW 


1.76 X 10" 


-7 


0.67 


39.00 


46.99% 




VO 


4.64 X 10" 


■10 


0.57 


38.28 


50.25% 




SMA 


4.94 X 10" 


-6 


2.33 X 10-^ 


50.94 


9.54% 


CI 


LMA 


4.70 X 10" 


-5 


0.38 


30.59 


82.99% 


Renorm. 


LOW 


1.99 X 10" 


-7 


0.77 


34.26 


68.57% 




VO 


4.61 X 10" 


-10 


0.59 


32.14 


77.36% 



Table 4: Best fits to the combined rates and SK day-night spectrum data in 
terms of oscillation into active neutrino with a fixed Xb — 0.75. 





Nature of 
Solution 


Aw? 
in eV^ 




tan^^ 


Amin 


g.o.f 




SMA 


5.28 X 10" 


-6 


3.75 X 10"^ 


48.39 


14.40% 


CI 


LMA 


4.65 X 10" 


-5 


0.49 


38.90 


47.44% 


Obsvd. 


LOW 


1.74 X 10 


-7 


0.71 


39.91 


42.95% 




VO 


4.55 X 10- 


10 


0.44 


37.17 


55.36% 




SMA 


8.49 X 10" 


-6 


1.78 X 10"^ 


50.77 


9.82% 


CI 


LMA 


4.64 X 10" 


-5 


0.51 


34.48 


67.61% 


enorm. 


LOW 


2.09 X 10" 


-7 


0.81 


33.47 


71.97% 




\'0 


1.59 X IQ- 


-10 


0.53 


30.63 


82.86% 
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Tabic 5: Best fits to the combined rates and SK day-night spectrum data in 



terms of oscillation into active neutrino with Xb free. 





Nature of 
Solution 




in cV^ 




tan^e 


• 


g.o.f 




SMA 


0.51 


5.25 X 10- 


-6 


3.44 X 10"^ 


46.83 


15.41% 


CI 


LMA 


1.18 


4.73 X 10" 


-5 


0.33 


32.32 


72.89% 


Obsvd. 


LOW 


0.88 


1.75 X 10" 


-7 


0.67 


38.75 


43.57% 




VO 


0.70 


4.55 X 10" 


■10 


0.44 


37.24 


50.44% 




SMA 


0.48 


4.66 X 10" 


-6 


2.32 X 10-4 


46.18 


17.01% 


CI 


LMA 


1.15 


4.71 X 10- 


-5 


0.36 


30.32 


80.80% 


Renorm. 


LOW 


0.83 


2.03 X 10- 


-7 


0.79 


33.18 


69.17% 




VO 


0.75 


4.63 X 10" 


-10 


0.55 


30.56 


79.92% 



Table 6: The values of the ratios Rsk / -^sno ^sno/ -^sno best-fit 
values for the LMA, LOW and energy independent solutions for the renor- 
mahsed CI and Xb—1-0 case. Also shown are the predictions for the maximal 
mixing (Pge — 0.5) solution and the experimental value of Rsk/^sno- 









pC'C / dNC 
^SNOl ^SNO 


tdES I dCC 
^SKl ^SNO 


Expt. value of 

r>ES 1 TfCC 
^SKI ^SNO 


LMA 


4.7 X 10-^ 


0.38 


0.30 


1.36 




LOW 


1.99 X 10-^ 


0.77 


0.45 


1.19 


1.33 ± 0.13 


energy-indep endent 




0.47 


0.56 


1.13 


maximal mixing 




1.0 


0.5 


1.15 
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Fig 1: The region of effectively energy independent solution in Am? — 
tan^ 9 parameter space where the solar neutrino survival probability 
agrees with eq. (3) to within 10% over the range of Ga to SK energies. 
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Fig. 3: The 90, 95 and 99 and 99.73% CL. allowed area from the global analysis of the total rates 
from CI (observed and 20% renormalised) , Ga,SK and SNO (CC) detectors and the 1258 days SK 
recoil electron spectrum at day and night, assuming MSW conversions to active neutrinos. 
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